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Abstract. We consider the transport of non-interacting electrons on two- and three-dimensional random 
Voronoi-Delaunay lattices. It was recently shown that these topologically disordered lattices feature strong 
disorder anticorrelations between the coordination numbers that qualitatively change the properties of 
continuous and first-order phase transitions. To determine whether or not these unusual features also 
influence Anderson localization, we study the electronic wave functions by multifractal analysis and finite- 
size scaling. We observe only localized states for all energies in the two-dimensional system. In three 
dimensions, we find two Anderson transitions between localized and extended states very close to the band 
edges. The critical exponent of the localization length is about 1.6. All these results agree with the usual 
orthogonal universality class. Additional generic energetic randomness introduced via random potentials 
does not lead to qualitative changes but allows us to obtain a phase diagram by varying the strength of 
these potentials. 

PACS. 7I.30.+h Metal-insulator transitions and other electronic transitions - 72.I5.Rn Localization effects 
(Anderson or weak localization) - 64.60.F- Equilibrium properties near critical points, critical exponents 

- 7I.55.Jv Disordered structures; amorphous and glassy solids - 73.20.Fz Weak or Anderson localization 

- 05.45.Df Fractals 


1 Introduction 

More than 5 decades ago, Anderson [T] showed that ran¬ 
dom disorder can localize a quantum particle in space. 
This phenomenon, now called Anderson localization, and 
the corresponding Anderson transitions between localized 
and metallic phases have since attracted lots of experi¬ 
mental and theoretical attention (see, e.g., Refs. [2l[3l[4] 
for reviews). Remarkably, important features of Anderson 
localization, such as the existence of a metallic phase and 
the qualitative characteristics of the Anderson transition, 
depend on the dimensionality and the symmetries of the 
system, in close analogy to conventional continuous phase 
transitions. 

This leads to a symmetry classification of Ander¬ 
son transitions. Originally, three universality classes were 
identified, based on the invariance of the Hamiltonian un¬ 
der time reversal and spin rotations. These classes (orthog¬ 
onal, unitary, and symplectic) correspond to the Wigner- 
Dyson classification of random matrices [KlIH] . For exam¬ 
ple, the orthogonal universality class contains systems 
that are invariant under both time reversal and spin ro¬ 
tation. Later, this classification scheme was extended by 
including additional symmetries (for a review, see, e.g., 
Evers and Mirlin [4]). Within each class, the properties 
of the Anderson transition are expected to be univer¬ 


sal. Specifically, all eigenstates in one-dimensional (ID) 
and two-dimensional (2D) systems in the orthogonal class 
are expected to be localized, implying the absence of an 
Anderson transitiorQ In contrast, the states in three- 
dimensional (3D) systems in the orthogonal class can un¬ 
dergo a transition from localized to delocalized as energy 
or disorder strength are varied; and this transition features 
universal critical exponents. 

These universal properties hold for uncorrelated ran¬ 
domness; spatial disorder correlations can lead to different 
behavior. Some short-range correlations have been shown 
to produce extended states at specific, isolated energies in 
an otherwise localized system, for instance in the so-called 
dimer model [9]. In contrast, long-range correlations can 
lead to the appearance of a true metallic phase (and thus 
an Anderson transition), even in ID. Some of these devel¬ 
opments are reviewed by Izrailev et al. uni. 

In recent years, topologically disordered systems (lat¬ 
tices with random connectivity) have attracted particular 
attention because phase transitions in such systems fea¬ 
ture surprising violations of the expected universal behav- 


^ Note that extended states (the so-called Azbel resonances 
mm can be found at isolated values of the energy. As their 
measure is zero and as they appear randomly in the energy 
spectrum, they can usually be neglected. 
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ior. This includes transitions in Ising and Potts magnets as 
well as the contact process, all defined on random Voronoi- 
Delaunay (VD) lattices |ll[ll2[[T3lll4[|15j . Barghathi and 
Vojta m solved this puzzle by showing that the 2D ran¬ 
dom VD lattice belongs to a broad class of random lat¬ 
tices whose disorder fluctuations feature strong anticorre¬ 
lations and therefore decay faster with increasing length 
scale than those of generic random systems. Such lattices 
are ubiquitous in 2D because the Euler equation for a 
2D graph imposes a topological constraint on the coordi¬ 
nation numbers; however, examples in higher dimensions 
exist as well. The suppressed disorder fluctuations lead to 
important modifications of the Harris m and Imry-Ma 
p!8lll9[l2Q[[^ criteria that govern the effects of disorder on 
continuous and first-order transitions, respectively. 

These results immediately pose the question of 
whether or not the unusual features of random VD lattices 
also modify universal properties of Anderson localization. 
Grimm et al. [22] studied the energy level distribution of 
a tight-binding model defined on a random tessellation 
[23l[^ similar to a VD lattice. They found level repulsion 
indicative of extended states in the metallic regime. How¬ 
ever, to the best of our knowledge, a systematic finite-size 
scaling (FSS) study that would permit the unambiguous 
identification of the metallic and localized phases (and the 
Anderson transition between them) has not yet been per¬ 
formed. 

In this paper, we therefore investigate noninteracting 
electrons on 2D and 3D random VD lattices. We perform 
FSS based on a multifractal analysis (MFA) of the elec¬ 
tronic wave functions. Our results can be summarized as 
follows. In 2D, we observe localized states for all ener¬ 
gies. In contrast, the 3D system features two Anderson 
transitions between localized and extended states close to 
the band edges. The critical exponent of the correlation 
length takes the value u ^ 1.6/m agreement with the stan¬ 
dard orthogonal universality class. This implies that the 
unusual coordination number anticorrelations of random 
VD lattices do not lead to qualitatively different behavior 
compared to the well-known Anderson model of localiza¬ 
tion on regular lattices. 

The rest of this paper is organized as follows. In Sec. 
we introduce the random VD lattices and define the 
Hamiltonian. We also discuss FSS as well as the MFA of 
the electronic states. Section [3] is devoted to the results of 
our simulations for both 2D and 3D systems. We conclude 
in Sec. HJ 


2 Model and methods 

2.1 Random Voronoi-Delaunay lattices 

The random VD lattice is a prototypical system with topo¬ 
logical (connectivity) disorder. It can be viewed as a sim¬ 
ple model for amorphous solids, foams or biological cell 
structures. The random VD lattice is defined as a set of 
lattice sites at random positions together with the bonds 
that connect nearest neighbor sites. These neighbors are 



Fig. 1. Left: 2D Voronoi-Delaunay construction for N = 100 
randomly positioned sites (dots) placed in a square of lin¬ 
ear size L = 10. The grayed outer margins show the peri¬ 
odic continuation of the original system (inner square). The 
Voronoi cells (polygons) are colored arbitrarily. The lines con¬ 
necting lattice sites show the Delaunay triangulation. Right: 
MFA scheme for the same setup. The system is partitioned 
into boxes of linear size 1. The smallest box size /min = 1 is 
chosen such that each box contains one lattice site on average. 


determined by the VD construction [25| as follows. The en¬ 
tire system area (or volume, in the 3D case) is subdivided 
into disjoint polygons (Voronoi cells) containing exactly 
one lattice site such that each cell contains all points in 
space that are closer to that lattice site than to any other 
site. The Voronoi diagram is the complete set of these 
Voronoi cells. Lattice sites whose Voronoi cells share an 
edge (or a face, in the 3D case) are considered neighbors, 
irrespective of their real-space distance. In 2D, the graph 
of all bonds connecting pairs of neighbors consists of tri¬ 
angles only. It is called the Delaunay triangulation. The 
corresponding construction in 3D leads to the Delaunay 
tetrahedralization, a lattice consisting of tetrahedra only. 
We emphasize that random VD lattices are not bipartite, 
i.e., they cannot be divided into disjoint sublattices A and 
B such that each bond connects an A site to a B site. A 
2D example of a VD construction can be seen in Fig. 

In the following, we consider independent, uniformly 
distributed random lattice sites of density unity contained 
in a square or cubic box of linear size L. We employ pe¬ 
riodic boundary conditions. More details of the algorithm 
that we use to perform the VD construction are given in 
the Appendix. 

In a random VD lattice, the coordination number 
(number of neighbors) fluctuates from site to site. This 
topological disorder is illustrated in Fig. The corre¬ 
sponding coordination number distribution can be seen in 
Fig-i The average coordination number in 2D is exactly 
= 6. This is a consequence of the Euler equation for 
a 2D graph consisting of triangles only. In 3D, the average 
coordination number is given by (/^) 3 d = 2 + (48/35)7r^ 
15.54 [26]. In both dimensions, the mean coordination 
numbers are higher than those of the regular square and 
cubic lattices often used in numerical localization studies. 
The standard deviations of the coordinations numbers are 

^ 1.33 in 2D and « 3.36 in 3D. This means that 
the disorder is moderately strong. 
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Fig. 2. 40 X 40 section of a 2D random VD lattice. The color 
of each Voronoi cell represents the coordination number of the 
included lattice site. 
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Fig. 3. Probability distribution of the coordination numbers 
for 2D (left) and 3D (right) random VD lattices. The data are 
averages over 500 systems of linear size L = 120 for 2D and 
L = 24 for 3D. Here, the maximal coordination numbers are 
15 and 35 for the 2D and the 3D system, respectively. 






on 2D (left) and 3D (right) random VD lattices for the case 
of purely topological disorder, VP = 0. The data are averages 
over 500 systems of linear size L = 120 for 2D and L = 24 for 
3D. The insets show magnifications of the right band edges. 


close to a particular energy value for larger system sizes, 
we use a sparse matrix algorithm based on the Jacobi- 
Davidson method m- Here, we treat system sizes up to 
L = 2000 for 2D and L = 140 for 3D. 

Figure shows an overview over the DOS resulting 
from these calculations for purely topological disorder 
{W = 0). As the VD lattices are not bipartite, the DOS 
is not symmetric with respect to the energy E = {). While 
the DOS near the low-energy band edge increases rapidly, 
there is a pronounced tail on the high-energy side, in par¬ 
ticular in 3D. Interestingly, this tail does not stretch much 
beyond E = (k). States in the far tail live on rare large 
clusters of sites with above average coordination numbers. 
The fact that there are almost no states with E > (k) re¬ 
flects the strong disorder anticorrelations of the random 
VD lattices that prevent the formation of large rare re¬ 
gions m- We also note that the DOS in the 3D band 
tail fluctuates significantly. We believe this is a finite-size 
effect. 


2.2 Anderson model 


We now consider the motion of noninteracting electrons on 
a random VD lattice. We describe it by means of a tight- 
binding model having one (Wannier) orbital per lattice 
site. The corresponding Hamiltonian 

H = j2\i){j\ + J2^i\i){i\ ( 1 ) 

(bi> * 


is analogous to the Anderson model of localization. The 
first term describes the hopping of electrons between near¬ 
est (Voronoi) neighbors. The hopping matrix element is 
constant and fixed at unity. This term contains the topo¬ 
logical disorder. The second term represents additional en¬ 
ergetic randomness. The Vi are random on-site potentials 
uniformly distributed in the interval [—VP/2, VP/2]. We are 
mostly interested in the case VP = 0 for which the disorder 
in the system is purely topological. However, for compari¬ 
son with the usual Anderson model of localization, we also 
consider nonzero random potential strength VP. 

To compute the densities of states (DOS), we directly 
diagonalize the secular matrices. The system sizes are 
1/ = 120 for 2D and L = 24 for 3D, limited by computer 
memory. To find the eigenstates of the Hamiltonian 0 


2.3 Finite-size scaling 

The scaling approach is a phenomenological description 
of the behavior close to the critical point of a continuous 
transition. We consider a general dimensionless measure 
A{E,L) characterizing the electronic states as function of 
energy E and a characteristic length, e.g., the system size 
L. Close to the critical energy Ec^ it fulfills the scaling 
form [281129] 


A{E,L) = E 




( 2 ) 


The description contains the relevant scaling variable 
associated with the relevant exponent v as well as the 
leading irrelevant scaling variable associated with the 
irrelevant exponent —y. The dependence of the scaling 
variables on E can be described in terms of the expansions 


V2r = cj + V and 

^ nl 

n=2 

mi 

f2, = i + y^^oj^ 

Z—^ /^ I 


n=l 


(3) 


( 4 ) 
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where uo = {E — E ^)/is a dimensionless measure of the 
distance from the critical energy. The scaling function 




Qr^L-^y 


n=0 


7 ,! 




(5) 


jii). For multifractal wave functions, Pq is expected to be¬ 
have as a power of the normalized box size A = //L, with 
the scaling exponent (mass exponent) 


Tg{^,l,L) = lim 


\nPg{$,l,L) 
In A 


(9) 


is expanded into a Taylor polynomial of the irrelevant vari¬ 
able. The coefficients are the regular single-variable scaling 
functions 


Fn 



'^an,kTk{^^vL-) . ( 6 ) 

fe =0 


They are expanded into Chebyshev polynomials T/c(') of 
the first kind and the kth order. This general expansion 
depends on the four expansion orders Ur, ni, rur, and m[ 
yielding (ur + 1)(ni + 1) + rur + mi + 1 + 0(ni) free param¬ 
eters. Their values are found by weighted nonlinear fits to 
the numerical data. The weights are formed by the recip¬ 
rocal variances of the data points. This formalism allows 
us to simultaneously determine E^^ and y. 

In order to get error estimates of these values we use a 
Monte-Carlo method [30]. It consists in building at least 
10^ synthetic data sets by adding noise to the original 
values. This noise is created by Gaussian random numbers 
with a standard deviation a equal to the individual error of 
each data point. By fitting A to these synthetic data sets 
we obtain distributions for the parameter values. These 
distributions are usually Gaussian. Large deviations, e.g., 
superpositions of multiple Gaussians, are interpreted as 
instabilities. The corresponding regressions are not well 
defined and will be neglected. (This could be avoided by 
increasing the accuracy of the data, changing the number 
of parameters, or varying the initial conditions.) We note 
that this method constitutes an error-propagation calcu¬ 
lation for the random errors only. Systematic errors are 
not detected completely, but the infiuence of different ex¬ 
pansion orders can be identified. 


2.4 Multifractal analysis 

Multifractal behavior is a feature of eigenstates at critical 
points m- The MFA is based on a standard box-counting 
algorithm. The boxes are squares or cubes for 2D and 3D 
systems, respectively. The d-dimensional system of size 
is partitioned into boxes of size (see Fig. [^. Using the 
probability 


. (7) 

rGbox b 

to find the electron in the bth box, the measure 

Pgi^,l,L) = J2^^li^,bL) ( 8 ) 

b 

is constructed from its qth moment. Depending on q, this 
measure is dominated by boxes with either large or small 


The singularity spectrum f{aq) is the Legendre transform 
of Tg, and it comprises the scaling exponents of the fractal 
dimensions of all moments. A parametric representation 
of the singularity spectrum can be obtained by calculating 
the singularity strength 


a 


Q 


^ = lim 
dq A^o In A 


and the fractal dimension 


fq = qaq -Tq= lim 


In A 


with 


jji =__ 

Pq{^,l,L) 


( 10 ) 


( 11 ) 


( 12 ) 


For statistical purposes, we utilize the ensemble average 
of these exponents IllEaES]. In particular, we use the en¬ 
semble averaged singularity strength 


with 


CXf] 


lim 

A^O 


Aq{l,L) 

In A 


„ _ (EbMbluMb).?. 


(13) 


(14) 


to perform the finite-size scaling analysis outlined in the 
last section. (•)^ denotes the average over different eigen¬ 
states. When determining the error of Ag, one has to take 
into account that the numerator and denominator are cor¬ 
related with each other m- 

Only integer box ratios A ^ — L/l are possible when 
partitioning the original system without overlap. A gen¬ 
eral ratio can be used by employing the periodic boundary 
conditions to fold back into the real system any protrud¬ 
ing box parts that arise for noninteger A“^ = Ljl [3^135] . 
To obtain a uniform sampling, we are then required to 
average over all possible box origins. 

To extrapolate to the thermodynamic limit A —0 
in Eq. (13), we perform a linear least-squares fit of the 
numerator versus the denominator at fixed L (box-size 
scaling). The standard deviations a of Aq(l) are utilized 
to determine the inaccuracy of the regression values. To 
approach the limit A ^ 0, sufficiently small values of I 
have to be taken into account. However, for the smallest 
/, the data deviate from the power-law behavior because 
corrections to scaling that stem from the irrelevant scal¬ 
ing variables become important. Their infiuence on the 
estimated exponent values can be reduced by optimiz¬ 
ing the fit boundaries. We use A G [0.1,0.3] for 2D and 
A G [0.3, 0.5] for 3D. 
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Fig. 6. Eigenstates for different energies (from left to right: -3.27, -2.89, 0.00, 5.80, and 6.31) of the same 2D VD lattice. Plots 
a) and e) show the wave functions of the lowest and highest eigenvalue, respectively. The amplitude of the individual site is 
given by the coloring of the Voronoi cells. The color white of the color scale corresponds to a homogeneously extended state, 
the color black corresponds to a state localized on a single site. 



distance min distance max distance 


Fig. 5. Probability distributions of the distances between 
neighboring sites: all (left), minimal (middle), and maximal 
(right) distances to neighboring sites. The data are averages 
over 500 realizations of 2D (upper panels) and 3D (lower pan¬ 
els) VD lattices with L = 120 and L = 24, respectively. 

In random VD lattices, the cells have varying sizes; 
the resolution of the MFA given by the smallest box size 
/min is therefore arbitrary. We use /min = 1, such that the 
average number of sites per box is 1 (see Fig. [^. Since 
empty boxes affect the MFA results, only box sizes larger 
than the maximum distance between neighboring sites are 
included in the MFA analysis. This guarantees that there 
are no empty boxes. In both dimensions, we use / > 4 in 
accordance with the distributions shown in the right panel 
of Fig. 


3 Results 

3.1 2D system 

Figure shows five representative eigenstates of the en¬ 
ergy spectrum of a 2D system with L = 120 and W = 0. 
The wave function of the lowest eigenvalue (Fig. HI) is 
exponentially localized. The maximum amplitude is con¬ 
centrated at the single Voronoi cell with highest coordi¬ 
nation number (here n = 12). States on the upper band 



Fig. 7. Left: Scaling of Aq(l) with fixed L = 1000 and different 
g at F = 6.2 and IT = 0. (Individual symbols are only shown 
for data points with A < 0.1.) All data with A > 0.1 follow the 
linear regressions (solid lines) closely. The dashed lines indicate 
the bounds of A G [0,1, 0.3] used for the fit. Right: Multifractal 
spectrum for different energies with moments —2<q<2 
staggered by 0.1. 


edge also show localized behavior very clearly (Fig. |^). 
However, the highest amplitude is orders of magnitude 
smaller and, correspondingly, the localization length is 
larger. More generally, states of smaller energy are influ¬ 
enced by local fluctuations of the coordination number. 
States of higher energy are driven by interferences and the 
probability amplitudes behave as non-integrable Chladni 
figures [36]. The number of antinodes increases with de¬ 
creasing eigenenergy. In summary, this is a visualization of 
localization in a classical and quantum mechanical manner 
for the states close to the lower and the upper band edge, 
respectively. Towards the band center the amplitude fluc¬ 
tuations increase and a classification by visual inspection 
is not easily possible anymore. 

We therefore turn to the MFA to quantitatively charac¬ 
terize the eigenstates as function of energy. Specifically, we 
consider the FSS behavior of aq for systems from L = 60 
up to 2000 with 1000 states per data point. Figureshows 
the box-size scaling of Aq{l^L) used to estimate aq for 
several q. Whereas values at higher A show the expected 
power-law dependence on A, results for small / deviate 
systematically due to finite-size effects. Therefore, they 
are neglected as discussed before. The multifractal spec¬ 
tra f{a) resulting from the analysis are also shown in Fig. 
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Fig. 8. Scaling behavior of aq=Q{L^E) in 2D VD lattices for 
energies close to lower and upper band edge. Each data point 
is an average over 1000 eigenstates close to the respective E. 
Error bars show la intervals. 


[^for representative energies close to the upper band edge. 
For the energy = 6.2, we obtain a parabola-like shape as 
expected. Towards the band center the parabola shrinks 
(see E = 5.7) and approaches the limiting case for com¬ 
pletely extended wave functions, namely the single point 
(2,2). In the opposite direction, close to the band edge, 
the spectrum transforms towards the extremely localized 
limit, the points (0, 0) for negative q and (oc, 2) for posi¬ 
tive q. 

The behavior of aq=o{E, L) (the position of the max¬ 
imum of the parabola) is shown in more detail in Fig. 
For energy values close to the band edge, <^^=0 generally 
increases with L. This means that these states are local¬ 
ized. The strength of the localization depends on E. The 
pronounced localization at the lower band edge seen in 
Fig-i is reproduced here. Strongly fluctuating wave func¬ 
tions near the band center have a value close io a = 2 
for all system sizes within their accuracy. The envelopes 
of corresponding wave functions have a very large local¬ 
ization length (larger than the system size). Importantly, 
the curves for different L do not show a common cross¬ 
ing point, implying that there is no Anderson transition. 
All states are localized. (The seeming crossings at very 
large aq=o{E^L) in Fig. 13 can be attributed to numerical 
artifacts.) 

We have also studied the effects of additional potential 
disorder W = 2 and IT = 4. It results in a broadening of 
the DOS and an enhancement of the localization. How¬ 
ever, qualitative changes compared to IT = 0 were not 
found. All states are localized. 


3.2 3D system 

The data analysis for the 3D VD lattice proceeds analo¬ 
gously. We use system sizes between L = 20 and L = 140. 
In contrast to the 2D system, we observe two Anderson 
transitions induced by purely topological disorder. These 
transitions are located close to the two energy band edges. 
The corresponding FSS behavior of aq=Q is visualized in 
Fig. [9|for the transition near the lower band edge and in 
Fig.j^for the upper edge. Thus, localized states exist only 
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Fig. 9. Scaling behavior of Og=o in 3D VD lattices as function 
of system size L and energy E close to the lower band edge. 
Each data point is an average over 2000 eigenstates of different 
lattices close to the respective E. Error bars show Icr intervals. 
Lines are result of the scaling approach with nr = 4, m = 1, 
mi = 0, and rrir = 3 (see Tab. 1 for more details). The solid 
and dashed gray lines mark the estimated transition point and 
the corresponding la error interval, respectively. Lower plots 
show the distributions of critical parameters obtained from 10^ 
synthetic data sets. 




near the band edges. The broad central area of the energy 
band encompasses extended states. This corresponds to 
the fact that the topological disorder is only moderately 
strong and results in localization behavior similar to that 
of a weakly disordered regular Anderson model. 

As shown in Fig. [3 the singularity strength aq fea¬ 
tures a smooth energy aependence for the transition near 
the lower band edge. Thus, the FSS approach is applica¬ 
ble. The results are presented in Tab. 1. It can be seen that 
the critical parameters Ec^iy and y are influenced by the 
irrelevant scaling variable. We compare different expan¬ 
sion orders to demonstrate the stability of the regression 
results. In particular, we neglect the data of small sys¬ 
tems L < 50 and use a FSS approach without irregular 
expansion (ni = 0). Such regressions show a higher robust¬ 
ness when changing initial conditions. Taking into account 
all data, we estimate the critical energy Ec = —5.122(3) 
and the critical exponent u = 1.62(4). The results of re¬ 
gressions without the irrelevant scaling variable deviate 
slightly {Ec = —5.121(1) and u = 1.58(4)) because of the 
neglected systematic shift of the intersections described 
by the irrelevant term. The quality of all fits is very close 
to unity. This indicates that the errors of the original data 
points were overestimated. 

For the transition near the upper band edge, the aq=o 
curves shown in Fig. [^are very noisy. This is caused by 
the low DOS (see right panel of Fig. and finite size ef¬ 
fects. The existence of a transition can still be inferred be¬ 
cause aq^o decreases with L for the smaller energies, while 
it increases with L for the largest energies. This indicates a 
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E 


Fig. 10. Behavior of aq in 3D VD lattices as function of system 
size L and energy E close to the upper band edge. Number of 
considered eigenstates differs for different data points and is 
reflected in the error bars (Icr intervals). 



Fig. 11. DOS for the 3D VD lattice with additional potential 
disorder. The data are averages over 500 systems of linear size 
L = 24. 


crossing and thus an Anderson transition. However, a clear 
transition point cannot be determined from the available 
data. In particular, the values are insufficient to perform 
a scaling analysis. (Also note that for small systems, gaps 
in the DOS appear. Therefore the number of considered 
eigenstates varies from 10 to 8000 between data points, 
leading to strong variations in the error bars.) A rough 
estimate of the critical energy is E^ = 16.31. 

We now turn to the effects of additional on-site dis¬ 
order, i.e., W 7 ^ 0. As in 2D, on-site disorder leads to a 
broadening of the DOS with increasing disorder strength 
W, this is demonstrated in Fig. El To find the localiza¬ 
tion phase diagram, we have performed a number of cal¬ 
culations with either fixed E or fixed W. System sizes 
L = 50, 60, and 70 were used with 10 data points for 
each L to determine the critical points. Figure shows 
the resulting phase diagram. The region of delocalized 
states is asymmetric. For positive energies, it extends to¬ 
wards higher disorder with a maximal disorder strength 
_ 47 5 at = 22. Figure [T^ shows the localization 
transition near the upper band edge for fixed W = 2b. 
The additional potential disorder smoothes the density of 
states and suppresses the artifacts in the <^^=0 data seen 
for the pure 3D VD lattice. The FSS approach is thus 
applicable. We obtain Ec = 19.69(1) and u = 1.55(14). 

The transition induced by W at fixed energy E = 10 
is studied in detail in order to obtain a more accurate esti¬ 
mate of u for systems with both topological and potential 



Fig. 12. Phase diagram of the 3D VD lattice. Localized and 
delocalized (red shaded area) states are separated by the mo¬ 
bility edge (thick red line). The inaccuracy of the measured 
transitions (red dots) is smaller than the symbol size. The thin 
gray and black lines are contours of the DOS between 0.0 and 
0.1 staggered by 0.005. 



Fig. 13. Scaling behavior (left) of Og=o in 3D VD lattices 
with additional potential disorder IT = 25, as function of sys¬ 
tem size L and energy E close to the upper band edge. Each 
data point is an average over 1000 eigenstates of different lat¬ 
tices close to the respective E. Error bars show la intervals. 
Lines are result of the scaling approach with rir = 4, ni = 0, 
TTii = 0, and TTir = 2. The estimated exponent is z/ = 1.55(14). 
The solid and dashed gray lines mark the estimated transition 
point Ec = 19.69(1) and the corresponding la error interval, 
respectively. Right plots show the distributions of critical pa¬ 
rameters obtained from 10^ synthetic data sets. 

disorder. The data analysis resembles the analysis of the 
transition tuned by for IT = 0, as described above. The 
details are summarized in Tab. 2. We observe a critical ex¬ 
ponent ly = 1.59(6) in agreement with the result for purely 
topological disorder. 


4 Conclusion 

To summarize, we have studied the effects of topological 
disorder on Anderson localization. To this end, we have in¬ 
vestigated the wave functions of noninteracting electrons 
on random Voronoi-Delaunay lattices by multifractal anal¬ 
ysis and finite-size scaling. In two dimensions, there is no 
Anderson transition as all states are localized, even in the 
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Table 1. Estimates of the critical parameters Ec, and y with their standard deviations for the transition near the lower energy 
band edge. The indices 1 and u denote the lower and upper bound of considered sizes, respectively. Q is the quality of ht. A^p and 
Ny are the number of points and degrees of freedom, respectively. aq=o was analyzed in the energy range E G [—5.17, —5.07] 
with 2000 eigenstates per data point. We set mi = 0 always. 


E 

G E 

V 

Gv 

V 

Gy 

Li 

Lu 

Tlr 

rii 

TTlr 


Np 

Np 

Q 

-5.1207 

0.0005 

1.580 

0.0423 

- 

- 

50 

100 

4 

0 

3 

45.7 

126 

117 

1.0000 

-5.1206 

0.0005 

1.577 

0.0404 

- 

- 

50 

100 

4 

0 

4 

45.5 

126 

116 

1.0000 

-5.1220 

0.0010 

1.626 

0.0449 

2.188 

0.444 

20 

100 

4 

1 

3 

51.6 

189 

174 

1.0000 

-5.1220 

0.0010 

1.623 

0.0494 

2.157 

0.447 

20 

100 

4 

1 

4 

51.6 

189 

173 

1.0000 


Table 2. Same as Tab. 1, but for the transition at fixed E = 10. aq=Q was analyzed in the disorder range W G [38.5, 43.5]. 


W aw Gj, y Gy L\ Lu ni rur Q 


40.4181 

0.0301 

1.591 

0.0294 

- 

- 

50 

100 

2 

0 

2 

43.0 

150 

144 

1.0000 

40.4224 

0.0318 

1.586 

0.0315 

- 

- 

50 

100 

2 

0 

3 

42.8 

150 

143 

1.0000 

40.3287 

0.0772 

1.585 

0.0554 

1.871 

0.539 

20 

100 

2 

1 

2 

51.7 

225 

215 

1.0000 

40.2947 

0.0827 

1.607 

0.0592 

1.671 

0.581 

20 

100 

2 

1 

3 

51.2 

225 

214 

1.0000 


absence of extra potential disorder. Adding random po¬ 
tentials further enhances the localization. In contrast, in 
three dimensions, the topological disorder of the Voronoi- 
Delaunay lattice induces two Anderson transitions close 
to the edges of the energy band, with localized states in 
the tails and extended states in the bulk of the band. If 
extra random potentials are added, the region of extended 
states first broadens with the broadening density of states, 
but then it shrinks and vanishes at some critical random 
potential strength. 

All these qualitative features agree with those of the 
usual Anderson model of localization. This means that 
the anticorrelations of the topological disorder [16] do not 
affect the universal properties of Anderson localization. 
This also holds for the critical behavior of the localization 
transition in three dimensions. The correlation exponents 
found in the present paper, viz., u = 1.62(4) for purely 
topological disorder and u = 1.59(6) for combined topo¬ 
logical and energetic disorder, agree within their errors 
with high-precision results for the usual Anderson model 
of localization imiET]. 

Why is Anderson localization not (qualitatively) af¬ 
fected by the disorder anticorrelations of the Voronoi- 
Delaunay lattice even though other continuous and first- 
order phase transitions are qualitatively changed and vio¬ 
late the usual Harris and Imry-Ma criteria? In the systems 
in which these violations have been found |11[I121[T3[I14[ 
US], the coordination number directly determines the lo¬ 
cal distance from the transition point because the effective 
coupling strength is simply the sum over the effects of all 
neighbors. Anticorrelations of the coordination numbers 
thus generate anticorrelated random-mass (or random-Tc) 
disorder. Anderson localization is more complex. In par¬ 
ticular, quantum interference effects are crucial, at least 
away from the band edges, and these effects are not cap¬ 
tured by the coordination number alone. Note, however, 
that the electronic states in the band tails, where the lo¬ 
calization is mostly classical, do seem to be influenced by 
the disorder anticorrelations (see discussion at the end of 


Sec. 2.2). Clearly, more work will be necessary to fully 


resolve the effects of topological disorder on Anderson lo¬ 
calization. 

This work was supported by the NSF under Grant Nos. 
DMR-1205803. We acknowledge useful discussions with H. 
Bar ghat hi. 


Appendix: Algorithm for creating random 
Voronoi-Delaunay lattices 

Computing the Voronoi diagram or the Delaunay triangu¬ 
lation of a given set of points (lattice sites) is a standard 
problem of computational geometry, and many different 
algorithms are discussed in text books and the research 
literature (see, e.g.. Ref. [38]). Our algorithm follows a 
suggestion by Tanemura et al. [39] and is based on the 
remarkable “empty circumcircle property” of a 2D De¬ 
launay triangulation. It states that every triangle formed 
by the bonds (edges) of the Delaunay triangulation has 
an empty circumcircle, i.e., a circumcircle that does not 
contain any other lattice sites. This is illustrated in Fig. 
[T^ Analogously, a 3D Delaunay tessellation features an 
“empty circumsphere property”: Each of the tetrahedra 
making up the tessellation has a circumsphere that does 
not contain any other lattice sites. 

Our algorithm considers the lattice sites one by one 
and finds the list of its (Voronoi) neighbors. This is done in 
a two-step process: (i) We first identify candidates for the 
neighbors based on their distance. All sites within a dis¬ 
tance Rc from the given site are included in the candidate 
list. For optimal performance, the cutoff radius Rc should 
be chosen as small as possible without missing neighbors, 
reasonable values depend on the structure of the set of lat¬ 
tice sites under consideration (see below), (ii) From these 
candidates, we then construct all triangles (in 2D) with 
empty circumcirles for which the given site is one of the 
vertices. In 3D, we construct all tetrahedra with empty 
circumspheres for which the given site is a vertex. The 
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Fig. 14. Left: Valid Delaunay triangulation, the circumcircle 
of the three chosen (red) sites does not contain any other sites. 
Right: This triangulation of the same set of sites violates the 
empty circumcircle property because the circumcircle of the 
chosen triangle (formed by the red sites) does contain other 
sites. 


entire algorithm can be coded efficiently in less than 400 
lines of Fortran 90 in 2D and less than 500 lines in 3D. 

For the current project, we have applied these algo¬ 
rithms to sets of L? (2D) or (3D) uniformly distributed 
random sites of density unity contained in a square or cu¬ 
bic box of linear size L. The sizes range up to L = 2000 in 
2D and L = 140 in 3D, and we use periodic boundary con¬ 
ditions. We found that the necessary values of the cutoff 
radius are quite small because the bond-length distri¬ 
bution of the random VD lattice drops off rapidly with 
increasing distance, see Fig. (in 2D, the tail is approxi¬ 
mately Gaussian). Empirically, we found that = 5.7 in 
2D and R^ = 3.7 in 3D are sufficient to find all neighbors 
in all the lattices we considered. 

For system sizes L > R^^ the computational effort of 
our algorithm scales approximately linearly with the num¬ 
ber of sites. To give an example of the performance, finding 
the Delaunay triangulation of 10^ sites in 2D takes about 
30 seconds on an Intel core 15-3570 CPU while 10^ sites 
in 3D take about 3 minutes. 
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